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ABSTRACT 

We study spiral galaxies using a simple global disk model as a mean for approxi- 
mate determination of mass profiles. Based on rotation curves and the amount of gas 
(Hi+He) , we find global surface mass densities consistent with the measurements and 
compare them with B-band surface brightness profiles. As a result we obtain mass- 
to-light ratio profiles. We give some arguments for why our approach is reliable and 
sometimes better than those assuming ad hoc the presence of massive non-baryonic 
dark matter halo. With this model, we study galaxies NGC 7793, NGC 1365, NGC 
6946 and UGC 6446. Based on THINGS- rotation curve we study also galaxy NGC 
4536 and compare the results with those we published elsewhere for the same galaxy. 
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1 INTRODUCTION 

Reliable determination of total mass of spiral galaxies is 
a crucial step in the search for non-baryonic dark matter. 
In particular, laboratory experiments aimed at detection of 
dark matter particles, could succeed only if there is sufficient 
amount of dark matter inside the disk of our Galaxy in the 
Sun's vicinity. 

Present models of mass distribution in spiral galaxies 
assume that in addition to stellar disk, bulge, luminous halo, 
there is massive halo of non-baryonic dark matter, most of- 
ten assumed spherical. It dominates the total galajcy mass 
and also governs the dynamics of outer parts of the stellar 
disk. This last feature makes it possible to test in some cases, 
if spherically symmetric rather than flattened dark matter 
distribution is actually present. 

Recently, we h ave analyzed mass distribution in sev- 
eral s piral galaxies (jjalocha et al.l [20071 . I2OO8I : iBratek et all 
|2008| ). We find that for some galaxies it cannot be spheri- 
cal at larger radii. Since we expect in this case a fiattened 
mass distribution to better approximate the gravitational 
potential at larger radii than the spherical one, we apply 
the global disk model. We then find that mass distribution 
of luminous matter accounts for rotation curves of the ex- 
amined galaxies. 

There is also increasing number of studies whose 
authors, using more involved galaxy models than so 
far, conclude, contrary to earlier findings, that lumi- 
nous matter accounts for rotation i n the i nternal galac- 
tic regions (ISellwood fc Kosowsk^ I2OO0I: lEvansI I2OOII: 



iPalunas fc Williamsl 2000l : Williams. Bureau fc Cappellaril 



120091 ). There is also possibility that rotation of galaxies in 
the outermost regi ons can be driven by magnetic fields, not 
by d ark matter (|Battaner. Garrido. Membrado fc Floridd 
119921 ). 

These results need to be seen properly in the context 
of non-baryonic dark matter searches. Currently, there is 
strong observational evidence of non-baryonic dark matter 
in clusters of galaxies, in particular, from gravitational lens- 
ing by clus ters and from coll iding clusters, as e.g. the Bul- 
let cluster (|Clowe et al.l l2006>). Also, the cosmological model 
provides support to the non-baryonic dark matter hypoth- 
esis, albeit at a more theoretical level. On the other hand, 
direct searches give n ull results, except for DAMA/LIBRA 
(|Bernabei et al.lT2007l ). however, these results are disputable 
and have not been confirmed by other searchers. 

The consistent picture might emerge if one assumed 
that clustering scale of non-baryonic dark matter is that of 
galaxy clusters. At smaller scales, non-baryonic dark mat- 
ter could be dilute. In particular, its contribution to spi- 
ral galaxy masses could be much less than anticipated in 
spherical halo models. It cannot also be excluded that spiral 
galaxies may differ from each other in abundance of dark 
matter - some of them could be reach in dark matter, the 
other might be devoid of it. 

Before presenting results for some spiral galaxies ob- 
tained with our model based on a global disk approxima- 
tion, we briefiy signalize some important issues concerning 
rotation curve modeling. 

The methods for ascertaining mass distribution in spi- 
ral galaxies commonly use very simplified, one would even 
say naive, one- or many-component galaxy models , that 
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assume symmetries of mass distribution (axial, spherical), 
properties of rotation fields (concentric circular), and com- 
position characteristics (constant M /L ratios, etc), which at 
most approximately agree with observations or sometimes 
even contradict them. Our model of spiral galaxies is also 
simplified, however, it possess some features which, accord- 
ing to us, gives it the advantage over other models - e.g. it 
does not assume in advance any M/L ratio or its constancy, 
but is obtained at output as a resulting M/L ratio profile, 
which seems more realistic. 

In the past, the prevailing opinion was that one needs 
dark matter inside galaxies to account for their rotation and 
to stabilize their allegedly cold disks, whereas there is in- 
creasing number of papers modeling the same galaxies with- 
out dark matter. This is achieved by considering sufficiently 
complex and thus also more "flexible" models, for example, 
those assuming wide, not thin, stellar disks, ect. There is 
also too much freedom in choosing galaxy models, and there 
is no general agreement which models are closer to reality. 

But, why to assume circular orbits, axial or spherical 
symmetry, constant mass-to-light ratios when this is not th e 
case for real galaxies? As is pointed out by Bosm^ (|l978h . 
most of spirals are not differentially rotating axisymmetric 
disks with only tiny disturbances due to spiral arms. There 
are large scale deviations from the symmetry, oft en with 
very large amplitudes or even with warped HI layers (jBosmal 
[l97i). 

The predictions of the simple galaxy models might 
further be modified if one preferred to reconstruct whole 
Doppler images rather than rotation curves, the more the 
latter are only some interpretations of Doppler images ob- 
tained not without additional assumptions. Going even fur- 
ther, what would happen if, instead of using the simplified 
models, one used a sort of templates for inferring matter 
distribution in real galaxies based on their Doppler images 
and the observed qualitative characteristics like the number 
of arms, the presence of a bar etcl Such templates, relat- 
ing mass distribution and its qualitative shape features to 
simulated Doppler images could be catalogued by elaborat- 
ing a large number of galaxy-like virialized stationary sys- 
tems from n-body simulations. Wouldn't this approach be 
the right step toward realistic modeling of mass distribution 
in galaxies? To what extent could this approach change our 
understanding of spiral galaxies and current opinions about 
their composition? 

1.1 Dark matter inside galaxies 

The question how much dar k matter is p re sent in galaxy's 
interiors remains still open jBosmalllQQg ). l|Sellwoodlll99i ') 
ISellwood &: Kosowskvl l|200Cf ) point out that the better are 
the data used for constraining the properties of galaxy halos 
and the better is the quality of simulations, there is more 
difficu l ties w ith dark matter models on galaxy scales. Also 
lEvanj l|200ll ) says that the case for CDM, having «a number 
of stubborn problems*, is weak judged on the data from 
galactic scales only. 

This is especially evident from the still growing num- 
ber of examples of morphologically various galaxies whose 
rotation curves can be accounted for without dark matter 
hypothesis by using non-canonical, more realistic approach 
to mass modeling. For example, even though still assum- 



ing the simplifying unrealistic axi-symmetry and that the 
radial mass profile should follow that o f lumi nosity with 
constant M/L ratio. [Palunas fc Willianii (|200Cl|) were able, 
under the strict maximum (thick) disk hypothesis, that is, 
not including a dark halo, to reproduce - with I-band M/L 
ratio of (2.4 ± 0.9) • /lys and consistently with stellar popu- 
lation synthesis models - the overall structure of the optical 
rotation curves for most of galaxies in their sample of 74 
spiral galaxies with various surface brightness profiles and 
rotation curve morphologies. They link the poorest fits with 
non-axisymmetric features, like bars or strong spiral arms, 
that influence the brightness profiles, the determination of 
inclination, major-axis position angle, etc. This situation, 
however, might be cured by considering more realistic mod- 
els. 

That the interiors of spiral galaxies cannot be 
dominated by dark matter is also suggested by dy- 
namic al arguments. As noted by iPalunas fc Williams! 

if the average projected surface mass density of 
dark matter halo was greater than that of the opti- 
cal disk, the common insta bilities would be suppressed 
(|Athanassoula. Bosma fc Papa ioannou 1987|) together with 
lopsided modes i Rix fc Zaritskvlll995i ) appearing in many 
disk galaxies. In addition, as they point out, there are mod- 
els that reproduce features of rotation curves within the op- 
tical disks on different scales and that the models ra rely re- 
quire dark matter halos ( Kalna js 1983 ; Kent 1986 ; BuchhornI 
119921 ). 

These interesting observations of absence of 
dark matter inside galajcies are best expressed by 
IWilliams. Bureau fc Cappellaril (|2009l). They speculate, 
having succeeded in fitting some galaxies without dark 
matter, whether Occam's Razor argument would not be 
applicable. Why to introduce something that is not needed? 
The authors also form less authoritative statements about 
dark matter presence in galaxies than others. What they 
merely say is that dark matter presence is only statistically 
significant. 

The presence of a dark matter component that dom- 
inates the outer mass distribution in spiral galaxies, 
seems evident from their optical rotation curves which 
do not fall off sufficiently fast close to their optical edge 

and from their extended 
. However, in the context 



(I Rubin. Thonnard fc FordI 



HI rotation curves (|Bosmalll97l 



19781) 



of what has been said above, one can also conjecture that 
the halo would disappear if the galaxy models were more 
realistic - allowing for changing mass-to-light ratios, assum- 
ing non-axisymmetric mass features (such as spiral arms and 
bars), if the assumed extinction models was altered eic. Such 
models should be also taking into account that the veloc- 
ity field in galaxies is much more complicated than circular 
and that the published rotation curves are obtained by re- 
ducing a two-dimensional Doppler image of a galaxy into a 
single one- dimensional curve (such an image by itself car- 
ries already a reduced, partial information about rotation), 
and that this reduction is performed not without additional 
assumptions concerning the interpretation, averaging meth- 
ods, etc, of the Doppler image (for example, an averaging 
method which produces rotation curve by calculating pairs 
[r, y {rv'^)/r\ scanning a Doppler image, where the averag- 
ing at a given r is carried out over a window of size Ar and 
weighted by, say, a function of local intensity of the image 
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within the window, is presumably good for circular orbits of 
test bodies in a spherically symmetric gravitational poten- 
tial, but not necessarily for a self-gravitating system which 
may be flattened). 

As concerns dark matter in outer galaxy regions, there 
is also possibility, never taken into account in modeling of 
flattened galajdes, that rotation of slightly ionized gas at 
large radii may be driven also by magnetic fields whose 
strength may be comparab le with that of gravitation of 
whole galaxy at these radii |Sanchez-Salcedo fc Reves-Ruij 



I2OO4I : iKutschera fc Jalochall2004 '). helping to lower (or even 
eliminate) the amount of dark matter needed in the out- 
skirts of galaxies where rotation curves are usually flat 
iBattaner. Garrido. Membrado fc Floridolll992h . 

It is true that realistic modeling of galaxy rotation is 
very challenging and might be very difficult to achieve, but, 
at the same time it is not clear whether the naive models 
used so far, lead to trustworthy results - has the opposite 
ever been rigorously proved? 

1.2 Dark halo and stability arguments 

In modeling of galactic optical discs, it is assumed they are 
highly flattened, axisymmetric and cold. Ther e is a common 
opinio n linked with stability argument of Ostri ker fc Peeblei 
1 19731 ) ■ that for such disks to be stable within the optical 
radius, one necessarily needs the galaxies to be immersed 
within massive dark matter haloes. 

This conclusion comes about as follows. As observed 
in numerical simulations, when random stellar motions are 
small compared to the streaming rotational motion in flat- 
tened galaxies, without such a halo a bar would form in one 
rotation period, and only afterward ax;ial symmetry would 
be regained. During this stage, approaching stationary state, 
the disk would heat up and could not remain cold (this would 
not agree with the assumption that disks are cold). There- 
fore, cold disks must be immersed in (usually spherical) ha- 
los of sufficiently large mass. 

There are various criteria known for deciding whether 
a virialized rotati n g fiat tened system could be stable. 
lOstriker fc Peebles! l|l973h conclude that a spherical halo 
with halo-to-disk mass ratio of 1 to 2.5 and initial value 
of total rotational kinetic energy to total gravitational en- 
ergy ratio of 0.14 ± 0.03 renders the disk stable (of course, 
the halo mass beyond the disk's outer radius has no ef- 
fect and should not be taken into account). Their conclu- 
sions are supported by many computer simulations and also 
by linear stability anal ysis carried out by Toomre for Mes- 
tel disk (|Toomrelll964l ). Note, however, that for nonlinear 
differential equations (and hydrodynamical description of 
galax;ies leads to such equations) it is not a priori obvious 
that linear approximation appropriately describes solutions, 
the more that the gravitational potential also evolves along 
with matter, thus also the insta bility modes and the asso- 
ciated spectrum of eigenvalues. ISellwoodI l| 19831 ) criticized 
the criterion, as well as various other criteria, to be of lit- 
tle use. Nevertheless, he points out, they all lea d to simi- 
lar d emands o n mass distribution. As shown in l|Sellwoodl 
Il983h . based on I Fall fc Efstathioul j 19801 ) models, one comes 
to the conclusion that, for stability, total spherical mass in- 
terior to the half-mass radius of the disc sho uld b e more 
than twice the disc mass in the same region. iHohll (| 19761 ) 



and ISellwoodI l|l980l ) find similarly, using their numerical 
models, that for stability of differentially rotating disks the 
ratio of spheroidal mass to that of the disk interior to the 
outer radius of the disk, should be appro ximately 2 , while 
for rigid disks the mass could be smaller l|Hoh]||l976l ). This 
spheroidal-to-disk mass criterion is much more handy to use 
than the Ostriker-Peebles-like criteria in modeling galax- 
ies, since mass functions of different galaxy components are 
much simpler to find. 

It should be noted, however, that the assumption of 
cold disks is based on the belief that the local observation 
of small dispersion of velocities in solar neighborhood in our 
own galaxy can be extrapolated onto whole Galaxy and even 
onto all other flattened galaxies. E ven if the assumption was 
true, there is an argument due to iKalnaid (| 19871 ) that the 
halo is not so efficient a stabilizer as a bulge alone could be. A 
massive halo would change rota tion field of the bar in a wa y 
not measured in real galaxies (jPalunas fc Williamd[2000h . 
As for modest halo, it would be difficult to decide abo ut its 
presence since its effect would be small l|Kalnaislll987l ). It is 
also known that random motions of stars of the disc near the 
center have stabilizing property as well as a bulge or a thick 
disk. Neglecting this effect could result in overestimating the 
bul ge+halo mass re quired for global stability. 

SellwoodI (|l983l ) surmises that most galax;ies might have 



only moderate fractions of mass in the bulge+halo. As 
concerns the internal galaxy regions comprising the opti- 
cal disk, the amount of da rk matter is very poorly known 
jPalunas fc Williamsll20()ol ). On the other hand, contrary to 
the common opinion that only sphe roidal halos ca n help to 
stabili ze cold galactic disks, bulg es (jKalnais 1987) or thick 
disks (|Athanas fc Sellwood| 1987 ) are n ow considered to be 
more efficient stabilizers "l|Sackettlll99^ ). What's more, the 
presence of dark matter in large amounts, dominating the 
disk mass, would simply suppress the instabilities required 
to form the observed structures such as bars and spiral arms 
(|Athanassoula. Bosma fc Papaioannoulll987l ). 

Because spherical halo mass outside the optical radius 
has no effect on disk stability, the facts mentioned above, 
revers e the disk stability argument by Ostrik er fc Peebleg 
(| 19731 ). frequently used as a support for the need of 
a massive dark matter h alo within the optical radius 
JPalunas fc Williams! I2OO0I) . Here, it is worth to add that 
ISellwood fc Evans! ( 200 II ) construct a model disk galaxy 
with an almost flat rotation curve which is stable without 
presence of dark matter. 



1.2.1 stability arguments and mass models 

In modeling of mass distribution in galaxies, provided one 
has indeed a model that properly describes a given galaxy, 
one in principle should check whether the found mass distri- 
bution is stable. If disks were indeed cold, then for stability 
the spherical mass component should be sufficiently large. 
It seems however, that for estimating of the mass distribu- 
tion in galaxies, the stability issue is not so important for 
astronomers, the more that the galaxy models used in this 
respect are very simple, one could say even naive. Majority 
of them assume idealized substructures such as disks and 
bulges with simple parametric profiles, and assume symme- 
tries that are not observed (spherical and axial). For exam- 
ple, the spherical mass to disk mass ratio in model used by 



4 



Ide Blok et ai] ((200^ (THINGS) is much less than 2, e.g., for 
galaxies NGC 4736, NGC 3031. Despite the fact that a dark 
matter halo was included, the ratio would suggest that the 
found mass distribution is instable as it does not meet the 
stability criteria mentioned previously, that is, the disk-to- 
halo mass is too large. 

There is a big 'but', however. All mass models for flat- 
tened galaxies are merely some idealizations that do not take 
into account different processes that make real galaxies sta- 
ble. The role of the models is rather to roughly estimate the 
gravitational potential from rotation curves assuming circu- 
lar orbits and axial symmetry, and so on, and then to derive 
the corresponding luass distribution which one believes to 
some degree approximate the actual amount of matter in 
real, and otherwise stable galaxies. For example, maximal 
disk models give "effective" surface density which should be 
treated rather as a substitution for different galaxy compo- 
nents both those flattened and spheroidal, described alto- 
gether by a single function of radius; only after a suitable 
de-projection one would infer about various galaxy compo- 
nents, thus it would be at least exaggeration to apply to such 
models the stability criteria which says that axisymmetric 
cool disks are instable. 

Of course, in real galaxies, masses are neither dis- 
tributed in infinitely thin disks nor in ideally spherically 
symmetric structures, orbits are not circular, rotation field 
has a very complicated structure, and mass distribution does 
not have axial symmetry - axial symmetry is clearly broken 
by bars and spiral structure. Thus it would be at least in- 
appropriate to qualify or disqualify different models on the 
basis of stability arguments even if they do not contain a 
halo of dark matter, since the models, by construction, do 
not describe real galaxies and ignore many subtleties. The 
models are only slightly more accurate than would be a di- 
mensional analysis, which states that the amount of electro- 
magnetic energy radiated by a galaxy should be given by 
L = aDcV^I (Gfi), where Vc is some characteristic velocity 
derived from Doppler image, Dc is its radial size estimated 
based on the spatial range of the luminous stuff, /j, is some 
mass-to-light ratio derived on population analysis ground, 
and a is some unknown dimensionless factor related to the 
particular geometry of mass distribution and of rotation 
field. All parameters apart from a can be estimated from 
measurements. The models only better approximate the pa- 
rameter a by giving its concrete model, more or less accu- 
rate. For example, the estimation should be better when for 
an elliptical galaxy one uses a spherical rather than a disk- 
like model, or, for a spiral galaxy, a disk model rather than 
spherical. The approximation should improve when one dis- 
tinguishes spherical and disk-like subcomponents, further, 
one can add a width to the disk, and so on. But it would 
be hard to treat these models seriously when dynamics and 
stability of a galaxy is concerned (especially when maximal 
disk model is used which treats on equal footing flat and 
spheroidal internal mass components). 

Without realistic modeling such difficult problems as 
structural stability, cannot be answered. By realistic model- 
ing we understand finding, for a given galaxy, such a spatial 
rotation field and mass distribution for an n-body system - 
a computer model of this galaxy - which would be able to re- 
construct at least the observed Doppler image of the galaxy. 
Reconstruction of a published rotation curve only, is not suf- 



ficient, because different groups obtain for the same galaxy 
different rotation curves depending on the details of their 
procedures applied for processing Doppler images. Such im- 
ages are also often nonsymmetric, unlike simple galaxy mod- 
els, and lot of information about galaxy kinematics they con- 
vey is missing by the fact that the images are by construction 
a sort of projection of a three dimensional rotational field 
onto the observation plane. 

1.3 Luminosity 

M/L ratio is important for studying the formation and evo- 
lution of galaxies, thus already from this standpoint it is 
rather strange to assume it is constant from the start, while 
it would be more desirable to determine it as a function of 
radius. Unfortunately, the theoretical M/L for an individual 
galaxy is not well known even when assumed constant. Its 
derived value depends strongly on the assumed mass func- 
tion for the stellar population, on the star formation history, 
on metallicity and depends also on extinction. It is also af- 
fected by the assumed internal extinction law, by the errors 
in luminosity measu rements, and by the distance dimming 
l|de Blok et al.ll200ll ). This ratio is generally assumed con- 
stant separately for the bulge and for the disk. Accordingly, 
one tacitly assumes that the contribution of bulge and that 
of disk to the rotation curve should be controlled only by 
the light distribution via a single number found by best fits 
as if brightness profiles were to trace the overall features of 
rotation curves. However, the assumption of constant AI/L 
ratio may not be correct, since it requires that extinction 
and popu lation gradients across the luminous disk should 
be small (jPalunas fc WiUiamsl bOOd ). Constant M/L ratio 
seems also improbable since it would imply that composi- 
tion of a galaxy is homogenous. 

Overestimating the importance of constant M/L may 
change considerably and in uncontrolled way the contribu- 
tion of bulge and that of disk to the overall rotation, the 
more that optical rotation curves are not featureless, while 
luminosity is usually fitted by simple laws (exponential, etc) . 
In particular, this may lead to errors in reconstruction of the 
relation between mass density and rotation of the disk com- 
ponent, since rotation of a disk is a nonlocal functional of 
the density profile, sensitive to the detailed local structure 
of the profile. This in turn may lead to wrong estimates 
of mass profiles of other mass components. To give only 
an example, many spirals simply do not have exponential 
disks or any other given by a simple analytical law, and 
assuming otherwise, results in large discrepancies in disk 
scale len gths published by d ifferent authors for the same 
galaxies IjPalunas fc Williamsi,2000) . Closely related to the 
non-uniqueness of the estimates of mass profiles is the so 
called disk-halo degeneracy - the uncertainty in the relative 
contributions of different galaxy components to the over- 
all rotation; for unambiguous disk-halo-bulge decomposition 
one requires additional constraints on mass distribution, in 
particular, such as the mass-to-light ratios or the assumed 
mass profiles and shapes of dark haloes. In finding such con- 
straints may also be help ful measurements of gravitational 
lensing in spiral galaxies l|Maller et al.ll200ol ). 

Usual rotation-curve fitting methods have too many free 
parameters and are not unique, require addit ional assump- 
tions, sometimes quite arbitrary. For example. Ide Blok et al.l 
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l|200ll ) present disk-halo decompositions with several differ- 
ent assumptions for stellar M/ L ratio for the same galaxy, 
or assume that only the minimum disk is present and that 
the observed rotation curve is due entirely to dark matter, 
or assume that stellar M/L ratio is zero (!) and take into 
account contribution of the atomic gas (H I and He). In con- 
trast, our approach does not have the arbitrariness problem, 
since we predict, not assume mass to light ratios. 

As pointed out bv lCiarduUo. Durrell. at "all (|2004l ) . there 
is only li mited evidence for constant M/L. What's more, the 
results of [Herrmann and CiarduUd (120051 1 for M33 suggest 
that this assumption is not valid at all, M/L ratio of the 
galaxy's disk increases approximately 5 times over 6 scale 
lengths of the inner disk. Their findings also suggest that 
there is very little dark matter in M33's central regions. 



2 MOTIVATION FOR THE USE OF A 
GLOBAL DISK MODEL 

As we have seen above, the case for dark matter inside galax- 
ies is very weak and the amount of ascertained mass of dark 
matter is strongly model-dependent. It should be stressed 
that the models are constructed in the same framework of 
Newtonian gravitation. There is increasing number of pa- 
pers where more and more complex models are considered 
compared with previous simple models of galaxies. They ac- 
count for the observed rotation of internal regions of flat- 
tened galaxies without dark matter, and some of them, based 
on maximal disk hypothesis, account for global rotation also 
without dark matter. The use of models with extended disk 
finds additional support from the fact that dark haloes may 
be very flattened rather than spherical. For flattened sys- 
tems, however, gravitation gets more complicated and new 
effects difficult to tackle with, such as influence of external 
masses on internal orbits, must be taken into account and 
which would be absent if external masses were distributed 
spherically symmetric. Also stellar mass-to-light ratio need 
not be constant. Assuming otherwise, galaxy mass models 
become very stiff leading, for example, to disk-halo conspir- 
acy. As we saw, this also leads to discrepancies between pre- 
dictions. 

As concerns the presence of dark matter in spiral galax- 
ies, we therefore prefer a more cautious approach than just 
to assume the presence of the (unobserved) dark matter halo 
from the beginning. We propose to admit dark matter, only 
if models that composed of baryonic matter distributed in a 
flattened disk (of stars and gas) and of more spherical bulge 
and stellar halo, fail to properly account for the dynamics 
of the disk measured through rotation curves. Such an ap- 
proach, we believe, could signiflcantly reduce the undesirable 
mod el-dependence of the amount of CDM in spiral galaxies 
(fjalocha et al. 2008; Bratek et al. 2008). 

Most models of galaxy formation predict that M/L 
should be a declining function of distance, e.g. the inner 
regions of a galaxy form first and contain older, higher M/L 
ratio star populations (jCiardullo. Durrell. at al l2004l ). As 
follows from the previous section it is not correct to assume 
that the ratio is constant throughout galaxies. This is the 
reason why we prefer to determine mass distribution in disk 
based on measurements other than brightness profiles first, 
and only then determine mass-to-light ratios as a result of 



comparison of the obtained mass profile with the measured 
luminosity profile. In several galaxies we have studied so far 
using our approach we observe that this ratio indeed de- 
creases while density approaches that of hydrogen and he- 
lium known from measurements. This also suggest that dark 
matter is not needed for some galax;ies even in the outermost 
regions. 

Conceptually simple, direct method for determining the 
gross mass distribution in galaxies, is rotation curve inver- 
sion. Maximal disk model is an example. Provided rotation 
curve is known globally, the corresponding mass distribu- 
tion can be found via single integral functional. The method 
is independent of any conventional assumptions concerning 
decomposition of a galaxy into various idealized subsystems 
with their unknown (usually correlated with each other) pa- 
rameters determinable from best fits. Unfortunately, direct 
rotation curve inversion is not unambiguous if mass distri- 
bution does not have special symmetries, e.g. spherical sym- 
metry. 

Various ambiguities pertinent to gravitation of flattened 
systems, easily hidden by th e use of a particular model for 
the dark mass (jSackettll 19961 ). may influence the actual im- 
portance of different galaxy components. This fact seems 
not much appreciated. In addition, there is also ambigu- 
ity in mass-to-light ratios within optical disks, since ex- 
tended rotation curves alone give no indication of wher e 
the luminous disk could end (iBahcall fc Casertaiic Il985h . 
Consequently, as noted in (jPalunas fc Williamsll200ol l. there 
are mass models which flt rotation curves within the op- 
tical disk both with and without significant dark compo - 
nent (|van Albada fc Sancisil [l98^ : iLake fc Feinswod ll989D . 
What's more, rotation of a galaxy considered dark matter 
dominated in o ne model can be d evoid of dark matter in 
another model l|jalocha et al.|[2008h . 

For a fiattened mass distribution, local density is a non- 
local functional of the rotation field. And vice versa, velocity 
on a given orbit is affected both by gravitation of masses ex- 
ternal and internal to that orbit. For that example, it would 
be impossible to reconstruct dynamical mass function from 
rotation of luminous matter moving in a potential well of un- 
seen external very flattened halo of dark matter. Such dark 
haloes considerably flattened toward the st ellar plane, re - 
sembling a disk more closely than a sphere (jSacket^Igg^, 
with axis rat i os be twee n 0.1 and 0.3 where sugges t ed by 
ISackett et all j 19941 ) and lOllinel (| 19961 ). iBratek et all (|200^ ) 
analyzed errors resulting from extrapolation of the rotation 
curve beyond the last measured point and gave a criterion 
for their estimation in global disk model. These errors are 
the fundamental obstacle in using disk models. The recon- 
structed mass density is viable only out to a distance of 
60% of the radius of the last measurement point. The less 
serious is sensitiveness of t his method to noise cr i ticized in 
(|Binnev fc Tremaind [l987l ') . However, as ISacketd (|l997h il- 
lustrates, this leads only to 10 — 20% uncertainty. Therefore 
one needs a method to minimize such errors. 

Once one accepts that very flattened haloes of unseen 
matter can exist then one may lose any predictability. For 
example, since dark haloes are not seen, one may assume 
anything about them, say, that outside the last measurement 
point there is present a flattened dark matter ring of a given 
arbitrarily large mass and given shape, encircling the galaxy. 
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To balance its influence on the luminous matter, one needs 
correspondingly large amount of mass inside galaocy. This in 
turn would made impossible even to determine the mass- 
to-light ratio of the luminous matter since the ratio would 
depend on the unknown mass of the flattened dark matter 
halo. 

There is a simple criterion for deciding, whether the 
mass distribution at larger radii in some galaxy is flattened 
rather than spherical. In the presence of spherical dark mat- 
ter halo, the non-spherical component of mass distribution 
should be negligible. Then, only the radial component of the 
gravitational force is important. In this case Keplerian mass 
function G~^rv^{r) defined for rotation curve Vc{r) must be 
a nondecreasing function of radius 



_d_ 

dr 



>0. 



(1) 



In deriving this inequality, we assumed circular orbits of 
matter in the vicinity of the galactic plane, which is cus- 
tomary in modeling of rotation curves. For spherical sys- 
tems, the Keplerian mass function is identical with the true 
mass function. Thus, if sphericity condition ((l} is not sat- 
isfied at larger radii, the true mass distribution cannot be 
spherical, and spherical CDM halo is thereby excluded (as 
far as the usual assumptions in galaxy modeling are em- 
ployed), or other forms of dark matter are required forming 
flattened structures. 

We can use the global disk model especially for galax- 
ies with rotation curves breaking the sphericity condition at 
larger radii. Note, that the distribution of internal masses is 
not much important for the overall gravitation at larger radii 
(contribution from higher gravitational multipoles of inter- 
nal masses decreases quickly with radius), therefore, the disk 
model can be used also for the inner parts, in spite of the 
fact that the center may be rather spheroidal since, after 
de-projection of a part of the disk density onto a spheroidal 
component, this does not change the central mass signifi- 
cantly - this is an approximation, a standard procedure in 
the maximal disk model. Since rotation curves of disks satis- 
fying the sphericity condition are also possible (|Bratek et al.l 
one can try to apply the disk model also to spiral 
galaxies of which rotation curves do not violate this condi- 
tion. 



2.1 The model 

The disk model assumes that whole matter in a galaxy is 
distributed in the vicinity of the galactic equatorial plane. 
The model also assumes that orbits of stars are circular and 
that dispersion of velocities is negligible. The other simpli- 
fying assumption is that whole system is axisymmetric. The 
relation between surface mass density a{r) and velocity is 
then given by 



AGrV.p. I / a(x)^^dx 



(2) 



r (x^ 



where K and E are ell i ptic f unctions of the first and sec- 
ond kind (|Bratek et al.|[200^ . This integral is understood 
in the principal value sense [V.p.), since both summands 
are divergent, thus the integral must be skilfully calculated. 

In practice, however, it suffices to assume that mat- 
ter is located close to the galactic plane only in the galaxy 
outskirts where the external orbits are stabilized by gravi- 
tational potential of the galactic interior (this contribution 
to gravitational potential is almost spherical at large radii), 
while in the vicinity of galactic interior the density can be 
treated as an effective density understood rather as a col- 
umn density of masses projected onto the galactic plane - 
a substitution for otherwise spherical and stable bulges etc 
in internal regions. Interpreted in this way, this approxi- 
mated cold disc model of an otherwise stable spiral galaxy 
also avoids, or at least weakens, the stability counterargu- 
ments often raised against global disk models (see also the 
discussion in section [L2)) . 

It is seen form equation ([2]) that for a fiattened mass 
distribution of which the disk model is a limiting approx- 
imation, rotational velocity of matter on a given orbit is 
dependent both of masses distributed interior and exterior 
to that orbit. The inverse r elation is a functional of whole ro- 
tation curve, and it reads jjalocha et ahllioOSl : iBratek et al.l 
.20083 



'^{r) = ^— V.I 



(3) 



In the literature of this subject, for some mysterious reason, 
there is always used an (equivalent) expression containing 
derivatives of velocities w hich are hardly measurabl e with 
satisfactory accuracy, e.g. jBinnev fc Tremaind[l987l ). 

The crucial point now is that a rotation curve is never 
known globally, thus the mass distribution even in the re- 
gion where rotation is known, cannot be determined from 
rotation only. This is the manifestation of "non-locality" of 
gravity of flattened systems - one cannot determine their 
mass functions from the observed rotation and, therefore, 
additional constraints on this mass distribution is required. 
This is particularly important in the regions where mea- 
surements of rotation end, while the analogous error is neg- 
ligible in the internal galaxy regions. This property is fatal 
for reconstruction of mass distribution inside a galaxy when 
outside the luminous part of this galaxy there is present an 
external ring-like halo composed of dark matter; if it were 
spherical there would be no problem at all. 

Having realized the unpleasant properties of flattened 
galaxies we devised in (jjaloc ha et al. 2008) a method which 
may help to minimize this uncertainty or remove it. One 
does not have to assume anything beyond the last measured 
point, but simply use the available data of mass distribution 
as a lower bound for the amount of matter in regions where 
rotation for some reasons could not be determined, but the 
amount of gas is still measurable. Our goal is to examine 
how good is our model in reproducing the mass distribution 
in a sample of spiral galaxies. 

The thin disk model's role is to find such a global (j{r), 
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that the corresponding velocity of rotation v{r) defined by 
integral ((2|, would best overlap with the observed rotation 
curve Vc{r). This task, however, is not unique, since for ev- 
ery R, there is an infinite set of functions a{r) that give rise 
to the same v{r) for r < R, even though the (j(r)'s may 
significantly differ from one another for r > R and have dif- 
ferent total masses. To minimize this non-uniqueness, we 
look for mass distribution in flattened galaxies by itera- 
tions. A particular r ealization of such iteration proposed in 
IJalocha et al.ll2008l ') assumes that mass distribution of gas 
in a given galaxy is known for radii greater than the range 
of the rotation curve. 

We bring to the rea der's attention the fact, that 
mass distribution found in jjalocha et al.|[2008l ) for galaxy 
NGC4736 using this method, is consistent with all measure- 
ments, and this was possible without dark matter. T his re- 
sult i s interesting, since another model of this galaxy l|Kentl 
[l983) predicts almost 70% abundance of dark matter and 
does not well explain its rotation at large radii. It is there- 
fore natural to examine whether this galaxy is exceptional 
or, maybe, if there are other spiral galaxies with lower abun- 
dance of dark matter. 



3 RESULTS 

Below, we present results for several galaxies obtained 
with the use of global disk model. We app lied the itera- 
tion m ethod analogous to that introduced bv lJalocha et al.] 
l|2008t l. We are aware this is not a representative sample of 
galaxies. However, we think, that even when rotation of only 
a single galaxy could be explained without dark matter, this 
would be interesting as such and worthy to be presented. 

Initially, we planned to examine only galaxies with ro- 
tation curves breaking the condition of spherical mass distri- 
bution, because such galaxies are presumably flattened and, 
therefore, natural candidates to which global disk model 
is applicable. However, in order to see whether the model 
produces viable results for other galaxies, we included also 
galaxies NGC 6946 and UGC 6446 for which we also found 
good and complete data. 

As concerns the error analysis, our method by construc- 
tion exactly reproduces rotation curves, thus also fits rota- 
tion measurements within error bars. This is the reason why 
we show rotation curves without error bars. 



3.1 NGC 4736 

Here, we once again consider galaxy NG C 4736, which ear- 
lier was studied in l|jalocha et al.l [20081 '). This time, how- 
ever, we use a nothe r, newer rotation curve published by 
Ide Blok et al.l (|2008| ). In fi gure [ll this new curve is com- 
pared with that obtained bv lSofug (| 19971 ) we used last time. 
Both the curves are shown assuming distance 4.7Mpc dif- 
ferent from the previous 5.1Mpc. Although both curves do 
not overlap we come to comparable conclusions concerning 
the amount of non-baryonic dark matter. NGC 7436 still re- 
mains a galaxy in which the halo is not necessarily required. 
Surface density obtained in the global disk model smoothly 
approaches that of gas in the outskirts of the galaxy, we ob- 
tain also low M/L ratios (1.28 for I-band, 1.44 for V-band, 
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Figure 1. Results obtained in the global disk model for galaxy 
NGC 4736 at D = 4.7Mpc. a) Rotation curve: THINGS mea- 



surements [solid circles j i de Blok et al.l [20081 ). the model [solid 
line], and compar ison with high resolution rotation curve taken 
from l lSofu6|[l997t) . b) model global surface mass density [solid 
line], surface mass density of HI+He [open circles]; and sur- 
face brightness: B-band ] solid squares], V -band [solid triangles], 
I-band [solid circles I i Munoz et al.|[l989l) : c) mass-to-light ratio 
profile without inclusion of HI+He. 



and 2.01 for B-band). In addition, the local M/L ratio de- 
creases with radius for the largest radii. Therefore, despite 
the new rotation curve d oes no t brea k sphericity condition 
so strongly unlike that bv lSofuj l| 19971 ). we still consider this 
gala:xy as containing no dark matter or that abundance of 
dark matter is at most very low. 
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THING, de Blok et al. 2008 
A HI, Carignan & Puche 1990 
o Ha, 36cm telescope, Dicaire et al. 2008 
□ Ha, 3.6m telescope, Dicaire et al. 2008 
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Figure 2. Upper: Co mparison of the NGC 7 793 rotation curves 
obtai ned i n THINGS Jde Blqk et al.ll200^ . llCarignan Puchej 
ll99Cl) and l lDicaire et al.ll2008l) . Symbols stand for measurements 
and lines for model fits. Middle: surface mass densities obtained in 
global disk model compared with surface mass densities of neutral 
hydrogen and helium and with surface brightness in the B-band. 
Bottom: Mass-to-light ratios derived based on mass densities ob- 
tained for difl^erent rotation curves. 



3.2 NGC 7793 

We ta ke r otation curves publis hed by a) ICarignan &: Puch^ 
l|l990l ). b) Ide Blok et al.' ('2008'') and c) three rotation curves 
published byiDicaire et al. (2008). The latter three curves 
were merged such to obtain a single rotation curve. In ef- 
fect we have three distinct rotation curves for the same 
gala:xy. As seen in figure [2] the curves are similar which 
shows rotation measurements for this galaxy are trustwor- 



thy. The corresponding masses are also close to each other, 
1-1.02x101°Mq. Surface mass density from the global disk 
model smoothly converges to that of gas at the outskirts of 
the galaxy for all three cases. We obtain also low values for 
M/L rat ios (2.93 — 3 in B-b and). We stress that the at- 
tempt of Ide Blok et al.l (|2008l ) to account for the rotation 
with spherical CDM halo, gives very poor results. By giving 
up the halo, we obtain surface mass density in agreement 
with the amount of gas, and low M/L ratios. Therefore, we 
think it is reasonable to assume that the galaxy does not 
have any dark matter halo. 



3.3 NGC 1365 

Now we consider two rotation curves (jSofud Il997l ') and 
(|Zanmar et al. 2008) the latter suitably smoothed out in a 
way to conform with data points and error bars. We have 
noticed that surface mass density corresponding to the ro- 
tation curve has approximately exponential falloff (see fig- 
ure [H]). Therefore, for a while one may assume that surface 
mass density is a superposition of two exponential profiles, 
in order to check whether such mass distribution accounts 
for the observed rotation. We find, this is indeed the case, 
however the mass density becomes larger than that of gas 
(although quickly falling off). On the other hand, the sur- 
face density obtained with our method converges smoothly 
to that of gas already in a place where galaxy is still bright 
in B-band. This can be interpreted as that apart from HI 
and He there must be present other luminous stuff. In gen- 
eral, this is very interesting, that for this galaxy already the 
doubly-expone ntial disk alone suffic es to account for the ro- 
tation curve of (|Zanmar et al.|[2008l ). Irrespectively of which 
rotation curve is considered, the obtained mass-to- light ratio 
is low (2.91 to 3.15 in B-band depending on which rotation 
curve and method is used). This ratio falls-off with increas- 
ing radius in every case. As so, galaxy NGC 1365 can be 
modeled without non-baryonic dark matter. 



3.4 NGC 6946 



W e take three rota. t ion cu rves published by 



Sanders! (|l996h . 



bv lBoomsma et al.' (2008) and lde Blok et al. (2008') (the lat- 
ter were smoothed out). As is seen, the published rotation 
data differ si gnificantly from each other, particularly rota- 
tion curve of Ide Blok et al.l (|2008l ) exceeds much from the 
other. By using global disk model we obtained surface den- 
sities which smoothly converges to that of gas. Mass-to- light 
ratio is low based on rotation by Sanders (3.49 in B-band), 
but it becomes twice greater for rotation curve of THINGS 
(6.15) M/L ratio for THINGS grows slightly with radius, for 
the other two curves it is approximately constant. Summing 
up, based on available rotation data it is difficult to uniquely 
decide about presence of non-baryonic dark matter in this 
galaxy. The result depends of which rotation curve is used. 

Using our modeling method we can explain all these 
rotation curves. It is worth of pointing out that the ob- 
tained surface densities smoothly converge to that of gas. 
However, based on THINGS rotation curve, we obtain M/L 
ratio much greater than for NGC 4736, NGC 1365, NGC 
7793, similarly, for the other two rotation curves. 
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Figure 3. Upper: Comparison of the NGC 1365 rotation curves 
obtained in dZanmar et al .120081) and (|So^3l993). Symbols stand 
for measurements and lines for model fits. Middle: surface mass 
densities obtained in global disk model compared with a 2nd or- 
der exponential decay fit and with surface mass densities of neu- 
tral hydrogen and helium and surface brightness in the B-band. 
Bottom: Mass-to-light ratios derived based on mass densities ob- 
tained different rotation curves. 



Figure 4. Upper: Comparison of t he NGC 6946 r otation curves 
obtained in llBoomsma et al.ll2008h . l lSander jllQOeh and THINGS 
jde Blok et al.ll2008h . Symbols stand for measurements and lines 
for model fits. Middle: surface mass densities obtained in global 
disk model compared with surface mass densities of neutral hy- 
drogen and helium and surface brightness in the B-band. Bottom: 
Mass-to-light ratios derived based on mass densities obtained dif- 
ferent rotation curves. 



3.5 UGC 6446 

The AI/L profile of galaxy UGC 6446 grows with radius 
attaining values much higher than for the other examined 
galaxies, even after subtraction of HI and He. Such a sin- 
gular behavior of the M/L profile is expected at the edge 
of luminous disk, where the luminosity per unit of surface 
tends to zero, A(r) — > 0, while surface mass density re- 
mains finite. Then, the local AI/L ratio becomes singular, 



a(r)/\(r) — ^ oo. After subtraction of mass contribution of 
gas, the increasing M/L ratio indicates, that the stars at 
the outskirts of luminous disk, contain possible admixture 
of died stars, "black" dwarfs, neutron stars, or stars of too 
low a mass for nuclear burning (so called brown dwarfs). It 
is important to note, that the mass density of such a died 
star component needs only to be residual, such that total 
mass density is infinitesimally bigger than the mass density 
of shining stars, a, + e for the M/L ratio significantly in- 
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Figure 5. Upper: Comparison of the NGC 6446 rotation curve 
obtained in dVerheijen &c Sancis il l200ll) . Symbols stand for mea- 
surements and lines for model fits. Middle: surface mass density 
obtained in global disk model compared with that of obtained in 
spherical mass model, with that of neutral hydrogen and helium 
and with surface brightness in the B-band. Bottom: Mass-to-light 
ratio derived based on mass density obtained different rotation 
curves. 



creases at the edge: (cr* + e)/\ e/\ when both cr*(r) — >• 
and A(r) — >■ 0. This very small mass density, e, of died stars 
does not show up anywhere except at the very edge where 
it is strongly enhanced by small A(r). Contribution of died 
stars to the mass of galaxy is thus negligible. 

The following thought experiment provides a strong ar- 
gument against the existence of a spherical halo around 
galaxy UGC 6446. Suppose for a while that this galaxy 
is dominated by such a halo. Then, the mass distribu- 



tion could be considered spherical with good accuracy, thus 
the corresponding mass function could be calculated from 
M{r) = ^"^"^i (and this can be done, since the rotation 
curve of this galaxy satisfies the sphericity condition). Then, 
the corresponding column mass density a^oi. can be found 
by projecting the volume density onto the galactic plane by 
using the formula 



rv (r) 
G ' 



Ccol 



+ ■ 



+ hm f f /'•'-<3,, 



27vR^R-2-p2 



(4) 



where R is the distance of the outermost point on the ro- 
tation curve. Note, that our formula is very advantageous 
from the practical point of view, contrary to the (equivalent) 
handbook formula, as it does not require derivatives of veloc- 
ities, which are measured with very poor accuracy. We com- 
pare the column density acoi. with total surface mass den- 
sity obtained from the disk model and with its component 
surface mass density of gas (HI+He). The result is shown 
in figure [5] We see that close to the "edge", acoi. is lower 
than the amount of gas! We stress that the two surface mass 
densities seen in figure [S] correspond to two extreme ideal- 
izations - in the first the galaxy is approximated by a thin 
disc, in the other, by a spherically symmetric mass distri- 
bution. The realistic mass distribution must be somewhere 
in between, but, as is seen, closer to that in disk model. In 
view of these observations, we classify UGC 6446 as a galaxy 
without massive spheroidal dark halo. 



4 SUMMARY 

We considered 5 galaxies, among which one was reexamined 
with a new rotation curve. Our aim was to answer the ques- 
tion if the galajcies could be satisfactorily modeled without 
non-baryonic dark matter. In this regard we check whether 
the following conditions are met: the found mass distribution 
should account for the measured rotation curve, it should 
approximately converge toward that observed for gas (HI + 
He) in the outer parts of galaxies (as this is the lower bound 
for the amount of luminous matter), the mass-to- light ratio 
should be low and its profile should be a declining function 
of radius close to the galaxy "edge". 

We find that these conditions are satisfied by our re- 
suhs for galaxies NGC 4736, NGC 7793 and NGC 1365 
for all rotation curves. It is reasonable to state that these 
are galaxies with flattened, disk-like mass distributions, for 
which the disk model approximation works quite well. The 
case of galaxy NGC 6946 is not so clear. There are significant 
discrepancies between published rotation curves, and conse- 
quently, in the obtained M/L ratios. Galaxy UGC 6446 sat- 
isfies some of our criteria without non-baryonic dark matter, 
however its M/L ratio is highest and, in addition, it grows 
dramatically toward gala:xy "edge". In another test we as- 
sumed that the galaxy is dominated by a spherical dark mat- 
ter halo and, using a spherical model, obtained the resulting 
bulk mass density. This density, accounting for the galaxy 
rotation, is in this case lower than that of gas, which can 
be considered as a strong argument for absence of spherical 
non-baryonic dark halo. This argument does not exclude the 



11 



presence of a non-spherical halo. However, the mass needed 
in the galaxy outskirts to account for the rotation in the 
disk model is practically that of gas. It would be therefore 
worthwhile to consider other than CDM causes of the M/ L 
ratio growth. 

We have also signalized the controversial issue, to what 
extent the very simplified galaxy models like the ones com- 
monly used by astronomers or like that used in the cur- 
rent paper, can be treated as models of realistic galaxies. 
In particular, whether stability arguments can be applied to 
favourize or reject a particular galaxy model, since the role of 
these models, which serve only as an idealized approximation 
incorporating many assumptions, is to estimate different 
global characteristics of galaxies like total mass, the amount 
of gas, global mass-to-light ratio, dark-to-luminous mass ra- 
tio etc, not to reconstruct exact spatial mass distribution, 
the more that many effects like the influence of global mag- 
netic fields in a galaxy on the rotation of gas in the galaxy 
outskirts, are normally not taken into account, and these 
fields may be important for accounting for flat rot ation 
curves (jBattaner. Garrido. Membrado k. Floriddll992l ). 

Global disk model is a natural approximation for flat- 
tened galaxies such as spiral ones, as far as it is not assumed 
in advance, sometimes by force, that each spiral galaxy must 
necessarily be immersed in an invisible halo of dark matter. 
Moreover, we are able to find such mass distributions which 
perfectly account for rotation curves, while in models as- 
suming massive halo these rotation curves often cannot be 
satisfactorily reconst ructed, eg. in the case of NGC 7793 in 
(|de Blok et al.|[2008l ). or the dark halo found is too weak to 
account for stability criteria (NGC 4736 in (de Bl ok et al.l 
[2OO8:)) (whereas, historically, the role of such haloes was not 
only to explain flat rotation, but also to stabilize disks). Al- 
though our model is very simple, similarly as other galaxy 
models customarily used by astronomers, it is better by that 
it does not assume constant mass-to-light ratios but it pre- 
dicts the ratio as a resultant mass-to-light profile. Our re- 
sults pertain only to the particular galaxies examined by us 
so far. From the fact that there is no need of non-baryonic 
dark matter in galaxies NGC 4736 or NGC 7793, one cannot 
conclude that one does not need it in other spiral galaxies. 
To claim otherwise, further studies are needed. Nonetheless, 
we find interesting the existence of spiral galaxies without 
non-baryonic dark matter or galaxies with small abundance 
of such matter. 
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